ENUMERATION OF PERMUTATIONS BY NUMBER OF CYCLIC 
OCCURRENCE OF PEAKS AND VALLEYS 
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f"^ ■ Abstract. In this paper, we focus on the enumeration of permutations by number of cyclic 

occurrence of peaks and valleys. We find several recurrence relations involving the number of 
O ' permutations with a prescribed number of cyclic peaks, cyclic valleys, fixed points and cycles. 

Several associated permutation statistics and the corresponding generating functions are also 
studied. In particular, we establish a connection between cyclic valleys and Pell numbers as 
well as cyclic peaks and alternating runs. 
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1. Introduction 

Let [n] = {1, 2, . . . , n}, and let & n denote the the set of permutations of [n]. A permutation 

it G & n can be written in one-line notation as the word n = 7r(l)7r(2) • • • 7r(n). Another way 

of writing the permutation is given by the standard cycle decomposition, where each cycle is 

■^J- ■ written with its smallest entry first and the cycles are written in increasing order of their smallest 

£H entry For example, the permutation -ir = 64713258 G ©8 has the standard cycle decomposition 

^O: (1,6, 2, 4) (3, 7, 5) (8). 

A permutation -k = 7r(l)7r(2) • • • n(n) G & n is alternating if 7r(l) > 7r(2) < 7r(3) > 7r(4) < • • • . 
£NJ ■ Similarly, ix is reverse alternating if 7r(l) < 7r(2) > 7r(3) < 7r(4) > • • • . It is well known pQ that 

the Euler numbers E n defined by 



v^ n x 

^ y^E n — - = tan x + sec x 

a" 



J-i . n U - 

n=0 



count alternating permutations in & n . The first few values of E n are 1,1,1, 2, 5, 16, 61, 272, .... 
The bijection it i— > ir c on & n defined by ir c (i) = n + 1 — 7r(i) shows that E n is also the num- 
ber of reverse alternating permutations in & n . The study of the Euler numbers is a topic in 
combinatorics (see [22]). For example, Elizalde and Deutsch [9] recently studied cycle up-down 
permutations. A cycle is said to be up-down if, when written in standard cycle form, say 
(&i, 62 5 03, ■ ■ •)' we have 61 < 62 > &3 < • • • • We say that it is a cycle up-down permutation if it is 
a product of up-down cycles. Elizalde and Deutsch [9j Proposition 2.1] found that the number 
of cycle up-down permutations of [n] is E n+ \ . 

There is a wealth of literature on peak statistics of permutations (see [3J [H [TTJ [12l EH [T71 [23] 
for instance). For example, Kitaev |11] found that there are 2 n_1 permutations of [n] without 
interior peaks. Let ir = 7r(l)7r(2) • • • 7r(n) G <5 n . An interior peak (resp. interior valley) in ir is an 
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index i £ {2, 3, . . . , n — 1} such that 7r(i — 1) < ir{i) > 7r(i + 1) (resp. 7r(i — 1) > ir(i) < 7r(i + 1)). 
Clearly, interior peaks and interior valleys are equidistributed on & n . Let pk (ir) (resp. val(7r)) 
denote the number of interior peaks (resp. the number of interior valleys) in tt. An left peak in 
tt is an index i £ [n — 1] such that ir(i — 1) < ir(i) > ir(i + 1), where we take tt(0) = 0. Let 
pk (it) denote the number of left peaks in tt. For example, the permutation tt = 64713258 £ ©8 
has pk (it) = 2, val (tt) = 3 and pk (-k) = 3. 
For n ^ 1, we define 

W n (q)= Y, <? PkW and W n (q)= Y q PkHn) - 

It is well known that the polynomials W n {q) satisfy the recurrence relation 

W n+l {q) = (nq-q + 2)W n (q) + 2q(l - q)W' n {q), 

with initial values W\(q) = 1, W2{q) = 2 and W${q) = 4 + 2g, and the polynomials W n (q) satisfy 
the recurrence relation 

W n+1 (q) = (nq + l)W n {q) + 2g(l - q)W n (q), 

with initial values W x {q) = 1, W 2 (q) = 1 + q and W 3 (q) = 1 + 5q (see p2J A008303, A008971]). 
The exponential generating functions of the polynomials W n (q) and W n (q) are respectively given 
as follows (see [12]): 
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W(q, z) = l + Y W n (q)— = n h( Z. f . u ^= S- (1) 

t-^ n\ y/1 — qcosh{z^l — q) — smn(z- v / 1 — q) 

An occurrence of a pattern r in a permutation it is defined as a subsequence in tt whose letters 
are in the same relative order as those in r. For example, the permutation tt = 64713258 € ©8 
has two occurrences of the pattern 1-2-3-4, namely the subsequences 1358 and 1258. In [2], 
Babson and Steingrimsson introduced generalized permutation patterns that allow the require- 
ment that two adjacent letters in a pattern must be adjacent in the permutation. Thus, an 
occurrence of an interior peak in a permutation is an occurrence of the pattern 132 or 231. 
Similarly, an occurrence of interior valley is an occurrence of the pattern 213 or 312. Recently, 
Parviainen [14|, 115] explored cyclic occurrence of patterns over 6 n via continued fractions. 

In this paper, we focus on the enumeration of permutations by number of cyclic occurrence 
of peaks and valleys. The paper is organised as follows. In Section [21 we collect some notation, 
definitions and results that will be needed in the rest of the paper. In Section^ we present several 
recurrence relations. In Section 0J we discuss two classes of triangular arrays. In Section [SJ we 
compute two exponential generating functions of generating functions of permutations by their 
numbers of cyclic peaks/ valleys, cycles and fixed points. In Section [6l we establish a connection 
between cyclic valleys and the famous Pell numbers. In Section [71 we establish a connection 
between cyclic peaks and alternating runs. 
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2. Notation, definitions and preliminaries 

In the following discussion we always write a permutation tt £ & n in standard cycle decom- 
position. Let 

n = (c\, 4 . . . c\ x ){cl , 4 . . . 4) . . . (cj, 4 . . . <), 

and let a be a generalised pattern. Following Parviainen |144I15|. the pattern a occurs cyclically 
in 7r if it occurs in the permutation 

W (tt) = c 1 , c 2 , . . . c ix , c l , c 2 , . . . Cj 2 , . . . c x , c 2 , . . . c ik , 

with the further restriction that d, and c{ are not adjacent, where 1 < j < A; — 1. For 
example, 132 does not occur in (1)(2, 4, 5)(3), but 13-2 occurs exactly twice. An entry c\n in 
the cycle (c^,^, ■ ■ ■ cf.) is called a cyclic peak (resp. cyclic valley) of tt if c J rn _ 1 < c° m > c^n+i 
(resp. c 1 m _ 1 > dm < c^ +1 ), where 2 ^ m ^ ij — 1 and 1 ^ j ^ k. For example, the permutation 
(1, 6, 2, 4) (3, 7, 5) (8) has the cyclic peaks 6 and 7 and the cyclic valley 2. 

Let cpk (tt) (resp. cval (tt)) denote the number of cyclic peaks (resp. the number of cyclic 
valleys) of tt. The number of fixed points of tt is fix(7r) = ff{l ^ i ^ n : 7r(f) = i}. A fixed- 
point-free permutation is called a derangement. Let P n denote the set of derangements of [n] . 
Denote by eye (tt) the number of cycles of tt. For n ^ 1, we introduce the following generating 
functions: 

P n (q,x,y) = Y^ g^W^W/*^); 

T66 n 

K n («,X,y)= ^ g cval Wa .c y cW yfixW . 

7T66„ 

M re (g) = Yl Q cpkM =J2 M ^Q k l 

D n (q)= £ 9 cpkW = E^«^ fc ; 
£«(?) = E <f valw = E^ fc - 



Let 



p(n, t, s, r) = #{tt e <5 n : cpk (tt) = t, eye (tt) = s, fix (tt) = r}, 

and let 

v(n, t, s, r) = #{tt £ & n : cval (tt) = t, eye (tt) = s, fix (tt) = r}. 

The Stirling numbers of the second kind {^} is the number of partitions of [n] into k blocks. 
Let S n (x) = Y2=\ { n k) xk - !t is wel1 known (see [HI A008277]) that 

which is equivalent to 

S n+ i(x) = xS n (x) + xS' n (x). 
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The associated Stirling numbers of second kind T(n, k) is the number of partitions of [n] into k 
blocks of size at least 2. Let T n {x) = Y,k>\ T K k ) xk ■ li is wel1 known (see [HI A008299]) that 

T(n + l,k) = kT(n, k) + nT(n - 1, k - 1), 

which is equivalent to 

T n+1 (x) = xT' n (x) + nxTn^x). 

Clearly, P n (0,x, 1) = S n (x) and P n (0,x,0) = T n (x). For example, take a permutation 

7r = (vr(«i), . . .)(7r(«2), •••)••• (n(is), ■ ■ •) 
counted by p(n, 0, s, r). Recall that 

p(n, 0, s, r) = #{7r G 6 n : cpk (n) = 0, eye (n) = s, fix (n) = r}. 
Erasing the parentheses, we get a partition of [n] with s blocks. Hence 

^p(n,0,s,r) = j 

We say that a permutation ir is called a circular permutation if eye (71") = 1. Denote by C n 
the set of circular permutations of [n]. Each n G C n +i can be written uniquely as a cycle of 
the form tt = (1, 01, ai-, ... , a n ). Let </?(7r) = 6162 ■ ■ ■ b n , where bi = Oj — 1 for 1 ^ z ^ n. The 
correspondence (p : C n +\ >— > & n is clearly a bijection. Using the bijection 99, it is clear that the 
coefficient of x in the polynomial P n+ i(q, x, y) (resp. V n +i(q,x,y)) is W n (q) (resp. W n (g 1 )). In 
the next section, we present recurrence relations for the polynomials P n (q,x,y) and V n (q,x,y). 

3. Recurrence relations 

Theorem 1. The numbers p(n,t,s,r) satisfy the recurrence relation 

p{n + 1, t, s, r) = (2t + s - r)p(n, t, s, r) + p(n, t, s — l,r - 1)+ 

(r + l)p(n, t,s,r + l) + (n + 2-2t- s)p(n, t-l,s,r). 

Proof. Let n be a fixed positive integer. Let <7j G G n +i be the permutation obtained from 
a G & n by inserting the entry n + 1 either to the left or to the right of a(i) if i G [n] or as a new 
cycle (n + 1) if i = n + 1. Then 

, x |cyc(cr) if t G [n], 

eye (tr») = < 

I eye (cr) + 1 if i = n + 1; 

and 

fix (<t) — 1 if i G [n] and a(i) = i, 

fix(crj) = < fix (a) if i G [n] and <r(i) 7^ z, 

fix(o-) + l ift = n + l. 
Recall that 

jp(n + 1, t, s, r) = #{<Tj G & n+ i : cpk (ctj) = t, eye (uj) = s, fix {ai) = r}. 

It is evident that cpk (<tj) = cpk (cr) or cpk (<Tj) = cpk (a) + 1. There are four cases to consider. 
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(a) If a has t cyclic peaks, s cycles and r fixed points, then we can put the entry n + 1 on 
either side of some cyclic peak or at the end of a cycle of length greater than or equal to 
2. This means we have 2t + s — r choices for the positions of n + 1. As we have p(n, s, t, r) 
choices for a, the first term of the recurrence relation is obtained. 

(b) If a has t cyclic peaks, s — 1 cycles and r — 1 fixed points, then we can insert the entry 
n + 1 at the end of a to form a new cycle (n + 1). As we have p(n, t, s — 1, r — 1) choices 
for a, the second term of the recurrence relation is obtained. 

(c) If a has t cyclic peaks, s cycles and r + 1 fixed points, then we can insert the entry n + 1 
to the right of a fixed points. This means we have r + 1 choices for the positions of n + 1. 
This case applies to the third term in the recurrence relation. 

(d) If a has £ — 1 cyclic peaks, s cycles and r fixed points, then we can insert the entry n + 1 
in one of the n — 2{t — 1) — s = n + 2 — It — s middle positions and the last term of the 
recurrence relation is obtained. 

This completes the proof. □ 

We can express Theorem [1] in terms of differential operators. 
Corollary 2. For n ^ 1, we have 

P n+ i(q,x,y) = {nq + xy)P n (q,x,y) +2q(l - q) \- x{± - q)- 



+ (1-1/) 



dq dx 

dP n (q,x,y) 



dy 
In particular, 

„ s s„ / s , sdP„(q,x,l) , s dP„(q,x,l) , , 

P n+1 (q, x, 1) = (nq + x)P n (q, x, 1) + 2q(l - q) n ™ ' > + x(l - q) "^ ' ' . (2) 

By Corollary [21 we can easily compute the first few polynomials P n (q,x,y): 

Pi(q,x,y) =xy, 

P2{q,x,y) =x + x 2 y 2 , 

P3(q, x, y) = (1 + q)x + 3x 2 y + x 3 y 3 , 

Pi(l, x, y) = (1 + 5g)x + (3 + Ay + Aqy)x 2 + 6x 3 y 2 + x 4 y 4 . 

Recall that M n {q) = P n (q, 1, 1). For 1 ^ n ^ 7, using ([2|), the coefficients of M n (q) can be 
arranged as follows with M n ^ in row n and column k: 

1 

2 

5 1 

15 9 

52 63 5 

203 416 101 

877 2741 1361 61 
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Let D Ut k(q,x) be the coefficient of y k in P n (q,x,y). Clearly, the polynomial P n (q,x,0) is the 
corresponding enumerative polynomial on D n . Note that 



D n , k (q,x) = Yl 9 cpkW x cycW 



7T66 n 

fix (7r)=fc 



x k y^ a cpk ^ X CyC ^ 
kj ^ 

<rev„_ k 



l x k P n _ k {q,x,Q). 
k 



dP n (q,x,y)\ 

- D ni i(q,x) = nxP n -i{q,x,Q). 



Thus 

f dP n 

V dy , y=0 
Therefore, we get the following result. 

Proposition 3. For n ^ 2, the polynomials P n (q,x,0) satisfy the following recurrence relation 
P n+1 (q,x,0) =nqP n (q,x,0) + 2q(l-q) dPn ^ Xj °K x(l-q) dPn ^ x Xj0) +nxP n ^(q,x,0), (3) 
with initial values Po(q, x, 0) = 1, P\(q, x, 0) = 0, P-2(q, x, 0) = x and Ps(q, x, 0) = (1 + q)x. 

Recall that D n {q) = P n (q, 1,0). For 1 ^ n ^ 7, using ([3]), the coefficients of D n (q) can be 
arranged as follows with D nk in row n and column k: 



1 

1 1 

4 5 

11 28 5 

41 153 71 

162 872 759 61 

In order to provide a recurrence relation for the numbers v(n,t,s,r), we first define an op- 
eration. Assume that n is a permutation in 6 n with k cycles C\, C2, ■ ■ ■ , C k , where Cj = 
(c[, c?2, ■ ■ ■ , cj.) and 1 ^ j ^ k. Let <p : & n — > & n be defined as follows: 

. <KC u c 2 ,...,c k ) = (<KC 1 ),<i>(C2),...,<Kc k )). 

• For every cycle Cj, we have 

MA,---M i ) = (M),M),---><i>(4 i ))- 

• Let {01, 02, . . . , Oj.} be the set of entries of the cycle Cj = (cj, Cg, . . . , cj.), and assume 
that ai < 02 < • • • < a«.. Then 4>(a m ) = a,;. + i_ m , where 1 ^ m ^ i«. 

For example, 0((1,6, 2,4)(3, 7, 5)(8)) = (6, 1,4, 2)(7,3, 5)(8). Following [ID], we call this op- 
eration a switching. Clearly, if ir has £ cyclic valleys, then <j)(ii) has £ cyclic peaks, and vice 
versa. 
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Theorem 4. The numbers v(n,t,s,r) satisfy the recurrence relation 
v(n + 1, t, s, r) = (2t + 2s — 2r)v(n, t, s, r) + v(n, t, s — 1, r — 1) 

+ (r + l)v(n, t, s, r + 1) + (n + 2 - 2t - 2s + r>(n, t-l,s, r). 

Proof. Let n be a fixed positive integer. Let cij 6 S n +i be the permutation obtained from 
a £ <5 n by inserting the entry n + 1 either to the left or to the right of a(i) if i € [n] or as a new 
cycle (n + 1) if z = n + 1. Recall that 

v(n + 1, i, s, r) = #{c>i G 6 n+ i : cval (<Tj) = i, eye (<7j) = s, fix (oi) = r}. 

It is evident that cval(<Tj) = cval (a) or cval (cj) = cval (a) + 1. There are four cases to 
consider. 

(a) If a has t cyclic valleys, s cycles and r fixed points, then (ft(o-) is a permutation with 
t cyclic peaks, s cycles and r fixed points. Consider the permutation (j)(a), we can 
appending n + 1 either at the beginning or at the end of a cycle of length greater than 
or equal to 2. We can also put the entry n + 1 on either side of some cyclic peak of <j>(cr). 
This means we have 2t + 2s — 2r choices for the positions of n + 1. As we have v(n, s, t, r) 
choices for a, the first term of the the recurrence relation is explained. 

(b) If a has t cyclic valleys, s — 1 cycles and r — 1 fixed points, then we can insert the entry 
n + 1 at the end of a to form a new cycle (n + 1). As we have v(n, t, s — 1, r — 1) choices 
for <t, the second term of the recurrence relation is obtained. 

(c) If a has t cyclic valleys, s cycles and r + 1 fixed points, then we can insert the entry 
n + 1 to the right of a fixed points. This means we have r + 1 choices for the positions 
of n + 1. This case applies to the third term of the recurrence relation. 

(d) If a has t—\ cyclic valleys, s cycles and r fixed points, then we can insert the entry n + 1 
in one of the n — s — {s — r) — 2(t — 1) =n + 2 — It — 2s + r middle positions. As we 
have v (n, t — 1, s, r) choices for a, the last term of the recurrence relation is obtained. 

This completes the proof. □ 

We can express Theorem U] in terms of differential operators. 
Corollary 5. For n ^ 1, we have 

Vn+i{q,x,y) = [nq + xy)V n + 2q{l - q) \- 2x(l -<?)- 



+ (1 - 2y + qy) 



dq dx 

dV n (q,x,y) 



dy 
By Corollary [5l we can easily compute the first few polynomials V n (q,x,y): 

Vi(q,x,y) = xy, 

V2{q,x,y) = x + x 2 y 2 , 

V 3 (q,x,y) = 2x + 3x 2 y + x 3 y 3 , 

Vi(q, x, y) = (4 + 2q)x + (3 + 8y)x 2 + 6x 3 y 2 + x 4 /, 

V 5 (q, x, y) = (8 + 16g)x + (20 + 20y + 10gy)x 2 + (15y + 20y 2 )x 3 + 10x 4 y 3 + rr 5 y 5 . 
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Recall that M n (q) = V n (q, 1, 1). For 1 ^ n ^ 7, the coefficients of M n (q) can be arranged as 
follows with M n) k in row n and column k: 

1 

2 

6 

22 2 

94 26 
460 244 16 
2532 2124 384 

Let D n ,k(q,x) be the coefficient of y k in V n (q,x,y). Clearly, D Ujk (q,x) = (l)x k V n - k (q,x,0). 

Then 

'dV n (q,x,y)' 



= D nl (q,x) = nxV n -i(q,x,0). 
dy Jy= 

Hence we get the following result. 

Proposition 6. For n ^ 2, the polynomials V n (q,x,0) satisfy the following recurrence relation 

V n+1 (q, x, 0) = nqV n (q, x, 0) +2g(l-g) dVn ^ °) +2x(l-q) ^^ ° } + nxV n - 1 (q, x, 0), (4) 

TOi/i initial values Vo(q,x,0) = 1, Vi(<7,a;,0) = 0, V2(g, x,0) = cc and V3(g, x, 0) = 2a:. 

Recall that D n (q) = V n (q, 1,0). For 1 ^ n ^ 7, using (jl|), the coefficients of D n (q) can be 
arranged as follows with -D n) fc in row n and column &: 



1 

2 

7 2 
28 16 
131 118 16 
690 892 272 

4. On combinations of polynomials and Euler numbers 

Recall that M n (q) = P„(g,l,l), M n (q) = V n (q,l,l), D n (q) = P n (q, 1,0) and D n (q) = 
V n (q, 1,0). It is easy to verify that deg M n (q) = deg D n (q) = [^J anddegM n (<?) = deg D n (q) = 
L§J - 1 for n > 2. 

We define 

R n {q) = M n {q 2 ) + qM n (q 2 ) for n ^ 3. 
Let Rn{q) = Yl l k>o^-n,kQ k ■ Then for n ^ 3, we have 



M fe-i if k is odd, 

r). — = — 

n.k 



R •- = * n ' 



M n k if k is even. 

n >2 
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The nth Bell number B n counts the number of partitions of [n] into non-empty blocks (see |18^ 
A000110] for details). Recall that M„o = #{vr £ & n ■ cpk(7r) = 0}. Take a permutation 

tt= (ir(h), .. .)(Tr(i 2 ), ...)■■■ (ir(ij), .. .) 

counted by M n $. Erasing the parentheses, we get a partition of [n] with j blocks. Hence 

R n ,o = M n>0 = B n . (5) 

Set R\{x) = 1 and B,2(x) = 2 + x. For 1 ^ n ^ 6, the coefficients of R n {o) can be arranged as 
follows with R n ^k in row n and column k: 

1 

2 1 

5 6 1 

15 22 9 2 

52 94 63 26 5 

203 460 416 244 101 16 

We define 

In(q) = D n (q 2 ) + qD n (q 2 ) for n > 2. 
Let /«(<?) = Sfc>o^«: fc 9 • Then for n ^ 2, we have 

if fc is odd, 




if A; is even. 
Recall that -D n ,o = t^I 71 " £ ^n : c pk (vr) = 0}. Take a permutation 

vr = (vr(pi), . . .)(vr(p 2 ), ■■■)••• OOfc); ■ ■ ■) 

counted by D n Q. When n ^ 2, erasing the parentheses, we get a partition of [n] into A; blocks 
of size at least 2. Therefore, -D n ,o = S/c>i T(n, k) for n ^ 2, where T(n, k) is the associated 
Stirling numbers of second kind. Set h(q) = 1. For 1 ^ n ^ 6, the coefficients of I n (l) can be 
arranged as follows with I n y. in row n and column k: 

1 

1 1 

1 2 1 

4 7 5 2 

11 28 28 16 5 

41 131 153 118 71 16 

Note that 

M 2 k+i,k = D 2 k+\,k = #{vr S S 2 k+i ■ cpk (tt) = k} 
and 

M 2 k+2,k = R>2k+2,k = #{v e 5 2 fc + 2 : cval (vr) = /c}. 
Then from [9j Proposition 2.1], we get the following result. 

Proposition 7. For n ^ 0, we have R n+ i, n = I n +\,n = E n . 
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Let ao, Oi,G2 • • • jO n be a sequence of nonnegative real numbers. The sequence is said to 
be log-concave if <2j_iaj + i ^ af for z = 1,2, ... ,n — 1. We end this section by proposing the 
following. 

Conjecture 8. For n ^ 1, the sequences of coefficients of the polynomials R n (x) and I n (x) are 
both log-concave. 



5. Generating functions 

We present in the present section the exponential generating functions of P n (q,x,y) and 
V n (q,x,y). Towards this end, we define 



z 
n! 



P(q,x,y,z) := VP n (g,x,y)— and V(q,x,y,z) := Y] V n (q,x,y)'- 

Instead of computing P(q, x, y, z) and V(q, x, y, z) directly, we first consider P(q, x, 0, z) and 
V(q, x, 0, z) since their recurrence relations involve less partial derivatives. 

Theorem 9. The generating functions P = P(q, x, 0, z) and V = V(q, x, 0, z) satisfy the follow- 
ing partial differential equations (PDEs) 

dP dP dP 

(l-qz)- + 2q(q-l)- + x(q-l)-=xzP, 

dV dV dV 



and whose solutions are 
P(q,x,0,z)=e-* 

V(q,x,0,z) =e~ xz l q 
respectively. 



\/g-l \ ( VQ + cos(zy/q - 1) + Vg - 1 sm(zy/q - 1) 
VQ + 1 ) \VQ-cos(z^q-l) - y/q- 1 sin(zy/q - I) 

Vq^l 

yjq - 1 cos(zy/q - 1) - sm(zy/q - 1] 



2Vq 



Proof. We first prove the assertions for V = V(q,x,0,z). Multiplying §%§ by z n /n\, followed by 
summing over n J2 1, and using the fact that V\(q, x, 0) = 0, we get 



Y J V n+1 (q,x,0)- = Y J 



n>\ 



v f n \,o ft ,9V n (q,x,0) dV n (q,x,0) 

nqV n (q,x,0) + 2q(l - q) \- 2x(l - q)- 



Oq 



Ox 



+ nxV n _ 1 (q,x,0) 



111 



whence the PDE for V: 



dV dV dV 

: i-qz)- + 2q(q-l)- + 2x(q-l)-=xzV. 



(6) 



In view of the so-called /3-extension [20] (see also [191 §7] for an instance of it), we may assume 
that V(q, x, 0, z) = v(q, z) x for some function v. Then 

dV x _ x dv dV x _ ± dv dV 

— = xv x l — , — = xv x l — , — =v x ]nv 

oz oz oq oq oz 
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so that ([6]) becomes 

dv dv 

(1 - qz)— + 2q(q - 1) — = zv - 2(q - l)«lnv. (7) 

Next, we let w = Inf. Then 

1 dv dw 1 dv dw 

v dz dz' v dq dq 

and ([7|) becomes the following linear PDE: 

dw dw 

(l-qz)— + 2q(q-l)— = z-2(q-l)w (8) 

with auxiliary system 

dz dq dw 



1-qz 2q(q-l) z-2(q-l)w' 
Solving the ordinary differential equation arising from the first equality, namely, 

dz z 1 

dq~ + 2(q-l) = 2q(q-lY 
we obtain the characteristic defined by 

c\ = tan" 1 y/q — 1 — zy/q — 1, 

where c\ is an arbitrary constant. On this characteristic, ([8]) becomes the following linear ODE: 

dw tan" 1 y/q-l-ci 
2q(q - 1)_ = ^== 2(g - 1)«, 



whose solution is 



w{q,z) = ^TS~ ln ( ) +~ 

q 2q V 1 J Q 

= -- + - iJ* 1 ' 1 ) | /( tan ' 1 v / g T7T -^v^ 7T ) 



q 2q \ q J q 

where C2 = /(ci) and / is a function to be determined. 

Since V(q, x, 0, 0) = , u(g,0) :r = 1 for all x, it follows that v(q,0) = 1 so that 



/(tan" 1 ^— i) = il n ^-A_Y 



Letting u = tan * \/q — 1, we then have q = 1 + tan 2 u = sec 2 u and hence 

/ SGC 11 \ 

f(u) = -In — = = -In sin it. 

2 \secr u — 1/ 

Since 

/(tan -1 v 7 ?- 1 - z^- 1) = - In sin(zy/q - 1 - z \J q - 1) 



Vv 7 ? - lcos(zV^ - 1) - sin(z^q - 1)/' 



we obtain that 



F(g, x, 0, z) = u(g, z) x = e™ x = e"**/« 



Vtf- Icos^v^ - 1) - sin(z v / g- 1) 
By imitating the above calculations, one can obtain the corresponding assertions for P = 
P(q, x, 0, z), whose details are left to the interested readers. □ 



12 



S.-M. MA AND C.-O. CHOW 



Corollary 10. The exponential generating functions P(q, x, y, z) and V(q, x, y, z) have the fol- 
lowing explicit expressions: 

' y/q - 1 \ / y/q + cos(zy/q - 1) + y/q - 1 sm(zy/q - 1) 



P(q,x,y,z) = e xz ^-^ 
V(q,x,y,z)=e xz(y - 1/q) 



y/q + 1 / V y/q - cos(zy/q - 1) - y/q— I sin(z v / g - 1} 

Vq^l 

y/q - lcos(zy/q- 1) - sin(zy/q - 1] 



2y/q 



respectively. 



Proof. Let a £ & n have k fixed points. There are (?) choices of fixed points. If r is the 
partial permutation obtained by deleting all the fixed points of a, then eye (a) = eye (r) + k 
and cval (<r) = cval (r) since only ^-cycles (£ ^ 3) of a contribute to cval (a). Thus, V n (q, x, y) = 
Yl=o (T) v n-k(q, x,0)(xy) k follows. Hence, 



E^ x '4 = EEur^^>°)^!? 



?2>0 



n^0k=0 



E^Ew^,o) 



v n—k 



fc>0 



n^fe 



(n-Jfe)! 



e xyz e -xz/q 



Vq^l 



y/q- lcos{zy/q - 1) - sin(z^/g - 1) 
thus proving the assertion for V(q,x,y, z). The corresponding assertion for P(q,x,y,z) follows 
from similar consideration. □ 

In Section [H we combinatorially prove that the nth Bell number B n is the constant term of 
M re (0), i.e., 

Af„(0) = M„ i0 = B n . 
We now present a generating function proof of this result. Since M n {q) = P n (q, 1, 1), we have 

^M n (0)- = limVP n ( g) l,l)- 



n>0 



n! g-»o 



lim 

<?->0 



(i;-n 



•v/9 — 1\ { y/q + cos(zy/q - 1) + v 7 ? - lsin(zy/q - 1) 



i 

2x/q 



.y/q + l) \y/q- cos(zy/q - 1) - ^/g - 1 sin(,zV<7 - 1) 
Denote the limit on the right by L. It is easy to see that L is of the indeterminate form 1°°. So, 
by l'Hopital's rule, we have 

1 L- lim 1 1 f^~ ^ (VQ + cos ( z V(1- 1) + \/<? - 1 sin(zV? - 1 



g->o 2 y^ \ ^ + 1/ \ y/q — cos(zy/q - 1) - y/q - 1 sin(z-v/<7 - 1) 
cosh z + sinh z — 1 = e z — 1. 



Consequently, 






n>0 



the right side being the exponential generating function of B n , thus proving M n (0) = i? n . 

Let i = \J— 1. Note that cosh(x) = cos(ix) and sinh(x) = — isin(ix). Combining ([1]) and 
Corollary [101 we get the following result. 
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Theorem 11. We have 

V(q, x, y, z) = e xz{ y- 1/q) W{q, z) ? . 

6. A Pell number identity 
The Pell numbers P n are defined by the recurrence relation 

P n = 2P n _!+P n _ 2 , n = 2,3,..., 

with initial values Pq = and P\ = 1. Other known facts about the Pell numbers include the 
Binet formula: for n = 0, 1, 2, . . ., 

(1 + V2)" - (1 - V2) n 
2^/2 

and the exponential generating function: 

°^ y n e (l+V2)y _ e (l-V2)y 

\ Pn ri = 271 ' 

71=0 

The Pell numbers can be expressed as (— l)-evaluation of V n (q, 1, 0), as the next theorem shows. 
Theorem 12. We have for n ^ 1, 

;C(-i) cvai(ff) = iVi. (9) 

Proof. By setting a; = 1 and q = — 1 in V(<7, x, 0, z), we see that 

V(-l, 1, 0, z) = e z (V^2cos(z^/^2) - sin(z V /r 2)) (-2)" 1 / 2 
_ e z [{V2 + l)e~^ z + (V2-l)e^ z ) 

z n 
= l + £P n _l — , 



since 



z — ' n! /n '^ n 

n>l ,yu n>0 



"' = (1+a/2)s _ p (l-v / 2)^ 



2^ JO 

e z [(V2 - l)e^ z + (V2 + l)e~^ z ] 



1. 



2V2 
Equating the coefficients of z n on both sides, ([9]) follows. D 
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7. Relationship to alternating runs 

Let 7r = 7r(l)7r(2) • • • ir(n) £ & n . We say that tt changes direction at position i if either 
n(i - 1) < 7r(z) > vr(i + 1), or 7r(z - 1) > 7r(z) < 7r(i + 1), where i e {2,3, ... ,n — 1}. We 
say that n has k alternating runs if there are k — 1 indices i such that 7r changes direction at 
these positions. Let R(n, k) denote the number of permutations in & n with k alternating runs. 
Andre pQ found that the numbers R(n, k) satisfy the following recurrence relation 

R(n, k) = kR(n -l,k) + 2R(n - 1, k - 1) + (re - k)R(n -l,k-2) (10) 

for re, fc > 1, where i?(l, 0) = 1 and R(l, k) = for k ^ 1. 

For re ^ 1, we define R n (q) = Sfc=i •ftfaj k)q k . Then by (flu]) , we obtain 

Rn+2(q) = q(nq + 2)R n+1 (q) +q(l- q 2 ) R' n+1 (q), (H) 

with initial value Ri(q) = 1. The first few terms of i? n (g)'s are given as follows: 

# 3 (<z) = 2g + 4g 2 , 

J R 4 (9) = 2g + 12g 2 + 10g 3 , 
i2g( g ) = 2<? + 28g 2 + 58q 3 + 32«j 4 . 

There is a large literature devoted to the numbers R(n, k) (see [18j A059427]). The reader is 
referred to [U [131 EI] for recent progress on this subject. In a series of papers [5J [61 [7j, Carlitz 
studied generating functions for the numbers R(n,k). In particular, Carlitz [5] proved that 



Instead of considering ^fc=o -^( n + -*-' k)q n ~ k as Carlitz did, we shall study 

R{q,z) = V^+i(g)-r. 



n^O 



It is clear that R(q, z) = H(^, qz). Hence 



</ 



,'q- 1\ ( \Jq 2 - l + qs\n(zyjq 2 - l) x 



9+1/ V 1 — qcos{z^Jq 2 — 1) 
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Let u = sec 1 q, i.e., seen = q. Then 

\J q 2 — 1 + q sm(zy/q 2 — 1) \ / sinu + sin^-y/q 2 — 1) 



1 — qcos(zy / q 2 — 1) / \cosu — cos(zy / q 2 — 1) 



COS 



sin 



2/^-zyV-l 



1 + cos(n — zy 7 ^ 2 — 1) 
1 — cos(n — z\/g 2 — 1) 
_ q + cos(z^/q 2 - 1) + yg 2 - lsin(zy / g 2 - 1) 
g — cos(zi/g 2 — 1) — y 7 ^ 2 — lsin(zy / g 2 — 1) 
where in the second, third and fourth equality, we have applied the following results: 

, ■ t fl> T\ o ■ fu + zy/q 2 -l\ (u-z^Jq 2 - 1\ 

sin it + sm(zy q z — 1) = 2 sin I - I cos I I , 

l n 7\ o • ( u + z^Jq 2 -l \ . ( u- z^Jq 2 -l \ 

cosu — cos(z\/q z — 1) = — 2 sin I I sin I , 



i n — T\ o i( u - z Vq 2 - 1\ r ro - i( u ~ z Vq 2 - 1 

cos(ti — z\J q l — 1) = 2 cos I * 1—1 = 1 — 2sm ' 

cos(zy / g 2 — 1) + \J q 2 — 1 sin(zy / g 2 — 1 



cos 



(t4 - z\/g 2 - 1) 



respectively. Thus, we have 



m z ) = ( q ~ X \ / / 9 + cos ( 2: V / g 2 - 1) + yV -lsin(zy / g 2 - 1) \ 
U + V Vg-COs(zy / ^ 2 ^T)-y / ?^Isin(zy^ 2 ^T)y' 
Alternatively, by imitating the proof of Theorem [91 it can be shown that R = R(q, z) satisfies 
the following linear partial differential equation 

... 9. oR , n . oR „ 

with initial condition R(q,0) = R\(q) = 1 and whose solution is (]12|) . 
Combining (|12p and Corollary 1101 we get the following result. 

Theorem 13. We have 

P(q,x,y,z) = e x <y^R(^,z)^. 
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